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ABSTRACT 

Context. Historically, stellar model atmospheres with scaled solar abundances (for all elements heavier then helium) have been widely 
used for analysis of the atmospheres of chemically peculiar (CP) stars. However, in reality, atmospheres of CP stars demonstrate a 
variety of abundances, not necessarily scaled to the solar composition. 

Aims. We study the effects of individual abundance patterns on the model atmospheres of CP stars. The main purpose is to conduct a 
systematic homogenous study to explore the abundance parameter space occupied by these stars. 

Methods. We calculated a grid of the model atmospheres of A and B stars (logg = 4.0) for different effective temperatures 
(Jeff = 8000, 9500, 11000, 13 000, 15 000, 20 000 K) and chemical compositions. We used the LLmodels code to compute model 
atmospheres with individual abundance patterns, varying the following elements: C, Mg, Si, Ca, Ti, Cr, Mn, Fe, Ni, Sr, Eu and He. 
We compared the computational results for these peculiar model atmospheres with those of reference model atmospheres of the solar 
chemical composition. 

Results. We present a homogeneous study of model atmosphere temperature structure, energy distribution, photometric indices in 
the uubyfj and Aa systems, hydrogen line profiles, and the abundance determination procedure as it applies to CP stars. In particular, 
we found that Si, Cr and Fe are the main elements to influence model atmospheres of CP stars, and thus to be considered in order to 
assess the adequacy of model atmospheres with scaled solar abundances in application to CP stars. We provide a theoretical expla- 
nation of the robust property of the Aa photometric system to recognize CP stars with peculiar Fe content. Also, the results of our 
numerical tests using model atmospheres with one or several elements overabundant (Si and Fe by + 1 dex, Cr by +2 dex) suggest that 
the uncertainty of abundance analysis in the atmospheres of CP stars using models with scaled abundances is less than ±0.25 dex. If 
the same homogeneous models are used for the abundance stratification analysis then we find that the uncertainty of the value of the 
vertical abundance gradient is within an 0.4 dex error bar. 

Conclusions. Model atmospheres with individual abundance patterns should be used in order to match the actual anomalies of CP 
stars and minimize analysis errors. 

Key words, stars: chemically peculiar - stars: atmospheres - stars: abundances 



• i— i , 1. Introduction ble peculiarities (usually enhanced) on low-dispersion spectro- 

J\. • . 1 . grams. For Am/Fm, these elements are Fe-peak and heavier ele- 

^ ■ A fraction of at least 25% ( |Schneideill993|) of the upper main se- ments Note that there ^ no gapfj in degree of peculiarity or in 

quence stars is known as spectroscopically peculiar stars. Their temp erature limits between these subgroups themselves or with 

spectra show anomalously strong (or weak) absorption lines respect tQ normal stars In fact> their domains smoothly over- 

of some chemical elements in comparison to those of normal lap each Qth ^ often providing temperature dependant spectra 

stars with the same fundamental parameters. These stars are of- anomalies (abundance effects); however this does not mean that 

ten called chemically peculiar (CP) stars, implying that unusual they m related in terms of the origin of their peculiarities, 
chemical composition (element abundances) rather than other 
causes is responsible for the spectrum anomalies. Consequently, 

this means that abundances in the atmospheres of CP stars are Subsequently, photometric data and high resolution spec- 
substantially different (enhanced or depleted) compared to the a - a analysis in the visual and ultraviolet regions confirmed pri- 
solar ones. mary abundance peculiarities and revealed new features of CP 
Effective temperatures of CP stars range from 6500 K to stars. As a result, modern conven tional classification is based on 
about 25 000 K, corresponding to spectral types from late F to the presence of a m agnetic field dWolfflll983tlBorra et alj|1983h 
very ear ly B. The widel y used simplified classification of the CP IShoriin et alj2002l) and leads to two main subgroups of CP stars, 
stars by iPrestonl d!974l) includes (in order of increasing effec- The sequence of magnetic CP stars (mCP) from the lowest tem- 
tive temperature): CP1 (Am/Fm stars), CP2 (Si, Cr-Sr-Eu clas- perature is Ap and Bp (Eu-Sr-Cr-Si), He-weak (Sr-Ti, Si) and 
sic magnetic A and B stars), CP3 (Hg-Mn, B-type stars) and CP4 He-rich B-type stars. The sequence of non-magnetic stars (in the 
(He-weak, B-type stars). The nomenclature specified in brackets same 7/ eff order) is Am, Ap Hg-Mn, and He-weak (P-Ga) Bp 
reflects which elements are responsible for the principal visi- stars. The field strength observed in magnetic stars ranges from 
a few hundred gauss to several tens of kilogauss. Other small 
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stars) usually are not considered within the framework of this 
classification. 

1.1. Observable abundances of CP stars 

Am stars. The classic definition for Am stars by IContil (1970) 
states that the Am class stars show apparent underabundant Ca 
(usually measured by the Can A 3933 (K) line) and/or an over- 
abundance of the Fe group and heavier elements. The Ca defi- 
ciency may or may not be accompanied by underabundant Sc. 
Typical abundance values for elements in comparison to those in 
the solar atmosphere are from -0.5 to -l.Odex for C, Ca, Mg 
and Sc, and from +0.5 up to +1 dex for the iron-peak (as Ti, Cr, 
Mn, Fe) and heavier elements. This includes the rare earth ele- 
ments (REE) which are enhanced by up to +1 dex as well. The 
range of effective temperatures for this class is 7000-10 000 K. 

Ap stars. The spectra of Ap stars show enhanced lines of Si (up 
to +2 dex), Sr and Cr (up to +3 dex). Many stars show an excess 
of the iron-peak elements from +1 to +2 dex. The REE elements 
(with an accent on Eu) are typically greatly increased and their 
abun dance values range from +3 up to +6 dex (see for exam- 
ple jBonsack & Wolfil(l980l IStrasser et ail 120011; iBagnulo et al l 
120031) . The C, N and O elements are ty pically underabundant by 
about -1 dex dRobv & Lambert) IT990). The deficiency of He is 
found to be more substantial for Ap than Am stars and spans 
from -1 to -2 dex. For Ap stars the abundances of the Fe-peak 
elements are not obviously correlated with the REE anomalies. 
Effective temperatures range from 6500 to 16 000 K. 

Hg-Mn stars. The g roup of Hg-Mn stars shows overadundantHg 
from +3 to + 6 dex (j jomaron et al.ll 19991) and Mn that can reach 
up to +3 dex dDolk et al.ll2003l) . Their spectra can exhibit over- 
abundance of some other elements (such as R Ga up to +2 dex, 
and others) and deficiency in Al, Zn, Co and others. The Hg- 
Mn stars are also peculiar in H e which is deficient by as much 
as -1.5 dex (iDworetskyl l2004t) . These stars occupy a range of 
effective temperatures from 10000 to 15 000 K. 

He-weak stars. The He-weak stars have much weaker helium 
lines than those in spectra of stars with the same photometric 
colors and hydrogen line profiles. The helium deficiency spans 
from about factor of 2 (0.3 dex) to 15 (~1.2dex). The group of 
He-weak stars consists of three subgroups: the P-Ga subgroup 
continues the HgMn stars towards higher temperatures; the Si 
and Sr-Ti subgroups are hotter analogs of Ap stars adjoining 
to the He-rich stars. In that sense, the He-weak group approx- 
imately continues the peculiar element abundances of its cooler 
and hotter analogs (for example strong Cr lines). Typical values 
of T eB range from 1 2 000 to 1 8 000 K. 

He-rich stars. Finally, the group of He-rich (or He-strong) stars, 
which are main sequence stars, show abnormally high helium 
abundance (A^e ~ A^h), while their effective temperatures are 
21 000-25 000 K. This group is outside the scope of this work 
due to its extreme nature and will not be considered in the cur- 
rent study. 

The general values for chemical composition given above 
turn out to be even more complicated in reality because "to those 
who really familiar with the spectra, some stars within a given 
class are as different from one an other as they are from stars 
in a different class" (ICowlevlll993L about Ap stars) and "large 
abundance differences between stars of similar temperatures and 
peculiarity classes may occur and do n ot correlate in any obvi- 
ous way with other physical properties" ( Bonsac k~& Worffll980l 



about CP stars). That is why we will use observed abundances 
in our study as suggested reference values only. 

1.2. Spectra analysis and model atmospheres 

CP stars, on the average, are much slower rotators 
(v e sin/ < 50km s~') than normal stars (except for the He- 
rich stars), which makes possible detailed line profile studies, 
providing better understanding of the CP stars phenomena. 
Spectra of magnetic stars show photometric, luminosity and 
line profiles variations (accompanied by the magnetic field vari- 
ability with the same period), whereas those of non-magnetic 
stars do not (althoug h some exceptions apply, see for example 
lAdelman et alj|2002l about Hg-Mn stars). 

The explanation for such variations is well developed, and is 
known as the rigid (oblique) rotator model (where oblique refers 
to the magnetic axis, which is inclined with respect to the rota- 
tion axis) with a non-uniform distribution of chemical elements 
over the stellar surface. Moreover, many magnetic and perhaps 
non-magnetic CP stars show signs of non-uniform distribution 
of the chemical elements w ith depth (abundance stratification) in 
their atmospheres (see e.g. lWade et alfc OOl; Sa vanov & H ubrig 
I2003I 

Here we want to draw attention to the well known fact that 
model atmospheres are required to determine stellar atmosphere 
abundances, which in turn are essential input parameters for stel- 
lar evolution theory (stellar structure modelling). For example, 
recen tly revised solar photospheric abundances ( Asplund et al. 
2005) due to development of a new highly sophisticated solar 
model atmosphere influenced th eoretical helioseismological pre- 
dictio ns in a major way (see e.g. lGuzik et al.l2005l ; lAn tia & Basu 
120051) . This illustrates why model atmospheres are a crucial part 
of the stellar astrophysics. 

In the case of CP stars, the surface non-uniformity, verti- 
cal stratification, presence of the magnetic field, and the ex- 
cess line blanketing make the spectral analysis quite complicated 
and uncertain. The usual analysis assuming a homogenous at- 
mosphere models with scaled abundances and w ithout magnetic 
fields (e.g., the famous classic atlas9 models bv lKuruc3 [l993a) 
may result in incorrect conclusions if such effects actually take 
place but are neglected in stellar atmosphere modelling. 

Despite this concern, classic model atmospheres are widely 
used for a variety of stellar astrophysics problems. They have nu- 
merous applications to CP stars: determination of fundamental 
parameters, abundance analysis, stellar magnetic field geometry 
research, detailed line profile study (including full treatment of 
the Zeeman effect), and reconstruction of stellar surface proper- 
ties by the Doppler Imaging (DI) technique (including all Stokes 
parameters). 

The problem is that while the correct spectrum analysis pro- 
cedure is an iterative process of comparison between observa- 
tional data, theoretical prediction, and the following corrections 
to the input modelling parameters, historically, it has developed 
that construction of model atmospheres and spectrum synthe- 
sis are almost co mpletely separated procedures. The study done 
bv lLesterl (1 19901) is one of the rare attempts to analyse the ef- 
fects of this separation. He found a small systematic error affect- 
ing the determination of stellar abundances if the line opacity 
background is neglected during the spectrum synthesis proce- 
dure, while the same opacity is fully considered during model 
atmosphere calculation. 

In fact, currently almost all spectral line profile analysis (in- 
cluding Stokes analysis) and abundance determinations are per- 
formed on the basis of classic model atmospheres (like atlas9 
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models), without taking into account magnetic field effects, pos- 
sible vertical abundance stratification, or even the specific chem- 
ical composition other than simply scaled to the solar abun- 
dance table. An initially picked model atmosphere is used for 
the whole analysis, usually not being recalculate d in accordance 
with r esults of this analysis. An indicative note bv lAdelman et al.l 
(1995) says that we need model atmospheres "that more closely 
match the actual abundance anomalies" to understand observ- 
able phenomena of CP stars. 

However, it should be noted that several authors have con- 
structed and used in dividualized/peculia r models in their re- 
search. For instance, iKupka et al.l d2004l) u sed a technique of 
model atmospheres calculation developed by Pis kunov & K upka 
(2001) (who applied the Opacity Distribution Function method 
to CP stars) to investiga te the famous depre ssion around 5200 A 
in spectra of CP stars. IShulvak et all (T2004) used a model atmo- 
sphere with individual and stratified abundances to fit the 5200 A 
depression, energy distributi on and hydrogen l ine pr ofile simul- 
taneously for a hot CP star. lLeBlanc & Moninl d2004l) made a big 
step forward calculating models with stratifie d abundances in a 
self-con sistent way. The most recent papers bv lKhan & Shu lyak 
(2006a,b) revealed effects of the magnetic line blanketing due 
to the Zeeman splitting and polarized radiative transfer on the 
model atmosphere structure and observable characteristics of 
mCP stars. 

Moreover, recent high resolution stellar spectra analy- 
ses dWade et al.ll2006t ILehmannet al.l l2006t IStiitz et al.l l2006t 



iFossati et al.ll2007t Tsilvester et al. 200 7} iNesvacil et al.ll2007l) ~as 
well as photometric analyses ( Stiit z & Paunzenl2006l) of CP stars 
done with help of the LLmodels code show that it is practical 
for model atmospheres with individual abundance patterns to be 
routinely utilized for spectral study. 

The main purpose of this paper is to study basic effects of the 
individual abundance patterns on the model atmosphere struc- 
ture, observable characteristics and abundance analysis of the 
non-magnetic CP stars. While it is obvious that all features of 
CP stars should be included in a whole cycle of the stellar spec- 
tra modelling, no systematic quantitative homogenous research 
has been attempted on this question before. 

This study continues our long-term investigation of model 
atmospheres of CP stars, highlighting some particular problems 
that observers encounter. Because of the complexity of the theo- 
retical phenomena and difficulties of working with real spectra, 
we focus on one of many problems associated with CP stars at 
a time. We provide indicative results that may help to estimate 
how uncertain specific aspects of spectral analysis can be. 



2. Calculation of model atmospheres 

Model atmospheres with individual abundance patterns calcu- 
lated in this study were constructed by means of the latest ver- 
sion of the LLmodels cod e (version 8.4) which uses direct treat- 
ment of the line opacity dShulvak et al J 120041) and provides us 
with a powerful tool to calculate model atmospheres of both CP 
and mCP stars within the LTE (local thermodynamical equilib- 
rium) approach. The code has been successfully applied for a 
number of studies (see references in the previous section). 



2.1. Technique 

To cover the stellar parameter space appropriate for CP stars, 
we calculated a grid of model atmospheres with effective tem- 
peratures r eff = 8000, 9500, 11 000, 13 000, 15 000, 20 000 K, a 



Table 1. Number of spectral lines (Mines), wavelength inter- 
vals (/l rang e) and corresponding numbers of spectral wavelength 
points (Appoints) for different effective temperatures (T e s) used 
for the opacity calculation. Lines were selected with a selec- 
tion criterion .1% using Kurucz's atlas9 model atmospheres 
dKuruczl fl993al) . ODF tables for [M/H] = +1 and enhanced 
chemical composition [M/H] = +2 to insure conservative pre- 
selection. 



T c fi, K 




a points 


Af lines 


8000 


500-50 000 


495 000 


2 693 251 


9500 


500-50 000 


495 000 


2760456 


11000 


500-30 000 


295 000 


2 870 658 


13 000 


100-30 000 


299 000 


3 267 218 


15 000 


100-30 000 


299 000 


3 342 524 


20000 


100-30 000 


299 000 


3 462 336 



single but typical value of the surface gravity log g — 4.0 and for 
a variety of abundance values. 

All models were calculated for the whole range of effective 
temperatures even if some of them are outside the actual temper- 
atures associated with a given class of CP stars (see Sect. II. Q . 
For example, models with overabundant Fe and Cr were calcu- 
lated for T e ff up to 20 000 K. This way we cover the whole pa- 
rameter space with model atmospheres to get a complete picture 
of changes and to analyse for any interesting behaviour. 

As usual, a standard reference point to specify a 
peculiar chemical composition was set to the solar 
abundances denoted hereafter as [M/H] = 0.0 (where 
[X/H] = [X] = log(JV x /Mi) - log(JV x /AWo and M means 
metals, i.e. any chemical element except H and He). Also, if 
the designation in square brackets is used then we omit the dex 
units. The list of the adopted peculiar abundance patterns is 
represented in the next section (Sect. |2~2| >. 

We used VALD dPiskunov et alJfl^lKupka et alJI 19991) as 
a source of spectral line data, including lines that originate from 
predicted levels. In total we extracted about 22 million lines; we 
then performed a preselection procedure to eliminate those that 
do not contribute significantly to the line opacity. 

For the line preselection procedure we u sed model atmo- 
spheres calculated by atlas9 dKurucz|[l993al) with correspond- 
ing fundamental parameters from our grid using the Opacity 
Distribution Function (ODF) tables for the largest valu e of the 
metall icity [M/H] = +1.0 available on the Kurucz CD dKuruczl 
1993b). At the same time we adopted the overabundant chem- 
ical composition [M/H] = +2.0 to insure selection of as many 
spectral lines as possible. The preselection was performed once 
for each set of models with the same effective temperature (but 
various abundance patterns). The line selection criterion, which 
is the relation between the line and continuum absorption co- 
efficients at the center of each line, was set to 0.1% (ten times 
less than we usually use). For the number of spectral lines pres- 
elected, and thus involved in computing the blanketing, refer to 
Tabled] 

To calculate model atmospheres of our grid, we used ei- 
ther atlas9 model atmospheres with metallicities [M/H] = +1.0, 
or model atmospheres already calculated within this study 
(whichever was more suitable) as an initial guess of the model 
structure to iterate from. Primarily to obtain faster convergence, 
we tried to use models with a close abundance pattern to con- 
tinue from them to more peculiar model atmospheres. Some 
models for T e g = 20 000 K required use of the grey approxima- 
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Table 2. The table represents different abundances adopted for model atmosphere calculations in this study. All abundances are 
given on a logarithmic scale with respect to the solar composition (i.e. dex units). Some elements and their respective values are 
accompanied by other element-value pairs indicated in parentheses. This designation means that models were calculated not only 
for abundance values of the element specified in its column but also for combinations of these elements with abundance variations 
of other elements specified in parentheses. If "all" is specified then all values listed in the respective column are combined with the 
main element value outside the parentheses. The star symbol (★) marks the reference solar composition model, the asterisk symbol 
(*) indicates non-typical abundance values for CP stars (see Sect. 13. 1 .41 , and the diamond symbol (o) points to the models with a 
decreased amount of He in order to normalize the sum of all abundances. 



[M/H] 


He 


C 


CNO Mg 


Si (Fe or Cr) 


Ca 


Ti 


Cr(Fe) 


Mn(Hg) 


Fe 


FeSiCr 


Ni 


Sr 


Eu 


Hg 


0.0* 


-0.417 


-1 


-1 -1 


+ 1 (all) 


-1 


+2 


+ 1 (all) 


+ 1 


+ 1 


+ 1 


+2 


+2 


+2 


+5 


+ 1.0 


-0.894 


-2 


+2* 


+2 (all) 


+2* 




+2 (all) 


+2 


+2 


+2 




+3 


+5 




+2.0o 


-1.894 






+3 (all)o 






+3 (all) 


+3 (+5) 















tion to make an initial model atmosphere structure, to avoid tem- 
perature correction freezing in the upper atmosphere. 

A logarithmic Rosseland optical depth scale logr loss was 
used as an independent variable of atmospheric depth spanning 
from +2 to -6.875 and subdivided into 72 layers. Convection 
and magnetic fields were neglected and a zero value for the mi- 
croturbulent velocity y tU rb was adopted to exclude their possible 
influence and focus on the abundance pattern effects. 

2.2. Adopted abundances 

Different chemical compositions were chosen in accordance 
with the typical observable abundance patterns found in spec- 
tra of CP stars (see Sect. II. Q . Consequently, we adopted vari- 
ous abundance values for the following elements: He, C, N, O, 
Mg, Si, Ca, Ti, Cr, Mn, Fe, Ni, Sr, Eu and Hg. Note that it is 
obviously almost impossible to perform a clear analysis on the 
basis of a large variety of abundance combinations. Thus, we 
chose to analyse the influence of each element separately, as a 
preliminary stage of the study, to understand which elements de- 
serve more attention. Besides typical amounts of the elements 
attributed to CP stars we considered several abundance values 
outside of those typical limits (see Sect. 13.1.41 for a description) 
as well as scaled to solar chemical compositions in order to ex- 
tend our study. 

The detailed data about all chemical elements considered and 
their respective values assumed in this study are presented in 
Table|2] Some model atmospheres were calculated for combina- 
tions of several chemical elements (which have been found the 
most influential elements). Such cases are specified both in the 
header of the table and in the table itself in the parentheses. 

For He-weak models we used three values of the deple- 
tion Afae/iVtotai = 0.001 (-1.894dex), 0.01 (-0.894dex), 0.03 
(-0.417 dex), while the solar amount is equal to 0.0783 
(0.0 dex). The hydrogen abundance was accordingly modified 
(increased) to normalize the sum of all abundances. The M/7V to tai 
value for the rest of the abundances stays the same as in the solar 
atmosphere (though in the log(M/Nn) scale they are decreased 
by 0.022-0.035 dex). 

Model atmospheres with scaled solar abundances were cal- 
culated for the values of [M/H] = +1.0, +2.0 and for the solar 
composition ([M/H] = 0.0). The latter models are called here- 
after reference model atmospheres because they are used to per- 
form the comparisons made in this study. 

The actual solar amount of each ele ment adopted i n this 
study was taken from the recent work by lAsplund et alj ([2005, 
Table 1) and converted into the AVMotai scale, where the fraction 
of hydrogen is 0.92 1 . 



Some model atmospheres required a decreased amount of He 
so the sum of all abundances would be normalized (see elements 
in Table [2] marked with a diamond (o) symbol). This includes 
model atmospheres with abundances scaled up by +2 dex (He is 
0.002) and model atmospheres with Si overabundant by +3 dex 
(He is 0.0483). For the rest of peculiar model atmospheres the 
sum of all abundances differs from unity by less than 1% (i.e. the 
sum of all elements £ Ni = 1.00 + 0.01), thus we left hydrogen 
and helium abundances for them unchanged. 

Finally, note that while for model atmospheres peculiar in 
some chemical elements the results are pretty obvious in ad- 
vance, nevertheless we adopted the peculiar values of those el- 
ements and performed the whole set of calculations in order to 
provide clear evidence for this point within the framework of the 
used data and techniques. 

3. Numerical Results 

In this section we present the numerical results of the modelling: 
the effects of the individual abundance patterns on the model at- 
mosphere temperature structure, energy distribution, photomet- 
ric indices, hydrogen line profiles, and on the abundance deter- 
mination procedure. 

The results are described with respect to the reference model 
atmospheres. This way we can clearly identify peculiar features 
provided by various abundance values of different chemical ele- 
ments. 

To find a generalized characteristic of which element pro- 
duces larger changes (is more influential) we use measurements 
of the average and maximum deviations in the whole parameter 
space of effective temperatures and abundances. 

While for a particular chemical element, usually the largest 
(or the smallest) abundance value produces the maximum effect, 
the comparison among different elements is not so straightfor- 
ward. The main problem, making the comparison more complex, 
is that different chemical elements may have different abundance 
ranges assigned to them. As a result, an element which is found 
to be more influential than another one may in fact produce less 
effect than the latter for some abundances or effective tempera- 
tures. 

It should be also remembered that there are some non-typical 
abundance values for some elements (marked by (*) in Table |2]i 
which are treated separately (Sect. 13 . 1 ,4b before drawing a con- 
clusion about the influence of these elements. 

3.1. Model atmosphere temperature structure 

In order to understand the influence of different abundances and 
their combinations on the model atmosphere structure, and thus 
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on the abundance determination through analysis of the line pro- 
files of particular absorption lines, we investigate the tempera- 
ture structure differences between model atmospheres with pe- 
culiar abundances and those with the solar composition. 

To represent such differences we find that the per cent scale 
(AT/T) provides a better view of the values of the changes than 
the absolute one. At the same time, for better understanding, 
it is convenient to know the absolute value rather than relative 
values. We note that a 1% change corresponds approximately 
to 50-90 K for T eS = 8000 K or 1 00-200 K for T eff = 20 000 K 
depending on the optical depth (from surface layers down to 
log Tress = 0). While this approximation is rather rough and de- 
pends on the effective temperature and particular shape of the 
temperature profile, we find that it is quite useful. 

We decided to divide all chemical elements into three groups 
according to the level of changes in the temperature structure 
they show in comparison to solar composition models with the 
same fundamental parameters (the reference models). A given 
element falls into a certain group if all the model atmospheres 
solely peculiar in this element demonstrate temperature structure 
changes which do not overstep the limits assigned to that group. 

3.1.1. Small (<1%) changes: He, C, N, O, Mg, Ca, Sr, Eu, Hg 

In this group we consider chemical elements those influence on 
the model temperature structure is not more than + 1 % in com- 
parison to the corresponding reference models with solar chem- 
ical composition. 

Models of the He-deficient stars demonstrate about 0.2% 
(or 20 K) temperature decreasing on the average and this cool- 
ing does not exceed 0.5% (or 60 K) for the highest value 
(20 000 K) of the effective temperature. For r eff < 1 3 000 K such 
changes are noticeable only from the optical depth log r loss = -2 
upwards the atmosphere. The temperature difference between 
models with the highest and the lowest He deficient values 
(A^He/Motai = 0.001 and 0.03) grows a bit with increasing T e g but 
does not become larger than 0.1% (or 10-15 K) for the highest 
value. H ere we note that analogou s results were obtained in the 
study by iGlagolevskii et alj ((2006, Sect. 3) where He-deficient 
model atmospheres were calculated by D. Shulyak with help of 
the LLmodels code in a similar way to the present study. We 
found that the variation of the mean molecular weight due to He 
deficiency is partly responsible for these small changes in the 
model atmosphere structure (see Sect. |3.2.TT >. 

For model atmospheres with the deficient element C, tem- 
perature changes do appear for high effective temperatures 
(Teff > 1 1 000 K), and show an approximately constant shift 
(plateau) with respect to the solar model, with the temperature 
about 1% (or 80-100 K) lower above logT ross = 0. For low ef- 
fective temperatures (T e g < 11 000 K) the maximum deviation 
is only as much as 0.3% (or 30 K) in some layers. The differ- 
ence in the temperature structure between models deficient in C 
by -2 or -1 dex is less than 0.2% (or 20 K), with more deficient 
models showing lower temperatures. 

The combined influence of the elements C, N and O each 
decreased by 1 dex results in a very small effect as well. The 
temperature structure falls in between the structures of the C- 
deficient models (-1 and -2 dex) and its deviation from the solar 
composition model grows from a vanishing value for the lowest 
effective temperature up to 0.2% (or ~15-20K) for the highest 
r eff = 20000K. 

Negligible changes of about 0.05% (or 5 K) in the range 
of the optical depths from logT ross = -3 downwards the atmo- 
sphere are demonstrated by the models underabundant by 1 dex 



in the elements Mg or Ca. For the upper atmosphere layers 
(log r loss < -3) such changes are larger and amount to at maxi- 
mum 0.1-1%. 

The model atmospheres with overabundant Sr (+2, +3 dex) 
and Eu (+2, +5 dex) differ from the solar composition models 
only by about 0.2%, which is clearly within the range of the 
computational errors. Actually, the temperature structure of the 
models overabundant only by +2 dex in Sr or Eu is practically 
indistinguishable from that of the solar composition models. The 
weak influence also can partly be explained by the lack of data 
about lines of the REE. 

Finally, an enormously increased amount of Hg (+5 dex) 
produces almost no changes in the model atmosphere temper- 
ature structure. In fact, it has the smallest influence among all 
elements considered above: specifically, it gives a 0.01% (or 1- 
2 K) change. 

3.1 .2. Setting the error bar threshold 

Despite the fact that some changes in the temperature structure 
can be confidently attributed to calculation errors while some are 
real changes, we have set the 1 % level as a threshold limit and 
consider all deviations below this limit as insignificant and un- 
reliable results. Though the same code was used throughout the 
whole work, we note that a variety of numerical problems (e.g. 
the transfer equation solution) have different valid methods to 
solve them, and once such methods are replaced with others, 
one would easily expect changes in the temperature structure 
of order 1%. The accuracy of the line list data is also limited. 
Moreover, it seems highly unlikely that the observations can pro- 
vide us with data accurate enough to perform any analysis with 
a precision higher than 1%. However, these concerns do not af- 
fect the generally high internal precision of the study itself, since 
all the models have been calculated the same way with the same 
basic physics, atomic data, modelling assumptions, etc. 

3.1.3. Moderate (1-3%) changes: Ti, Mn, Ni 

Model atmospheres and corresponding chemical elements as- 
signed to this group produce changes in the temperature struc- 
ture that clearly exceed the 1 % limit for some effective temper- 
atures but are almost always below the 3% limit. The elements 
Ti, Mn and Ni, which belong to this group, are Fe-peak elements 
and can be overabundant in stellar atmospheres (see Sect. II . lb . 

For models with excess Ti and Ni, the main line forming re- 
gion does not show changes that are actually larger than 0.3% 
for Ti, and 0.35-1.2% for Ni. Precisely, the region with less 
than ~1% changes spans from the deepest layers up to about 
log t 10SS = -3 optical depth. Changes of more than 1 % occur 
only above that level (exceeding 3% for Ni in some layers). We 
found that the maximum changes in the temperature structure 
appear at 1 1 000 K for Ti, and at 20 000 K effective temperature 
for Ni. 

The element Mn produces more profound effects on the 
model structure, increasing the temperature in the line forming 
region by about 1-2% and steadily decreasing it by ~2% in the 
outer layers (logr ross < -3). This temperature behaviour is very 
close to that demonstrated by models with enhanced Cr and Fe 
(see Sect. l3.1.5l >. where two distinct regions of heating and cool- 
ing determine the overall temperature profile. In particular, Mn 
produces changes which are closer in shape to those produced 
by Cr (see Fig. [2] as an illustration) than those by Fe, although 
all of them generally are quite similar. 



6 



S. A. Khan and D. V. Shulyak: Theoretical analysis of the atmospheres of CP stars 





















T „ = 8000 K 
'••••«......„„....•<• eff 


-6 -5 -4 -3-2-1 1 2 



N. ■••.;........•••"".,....„...••••■'" \ 

\ • 














T = 13000 K 

eff 





-6 -5 -4 -3 -2 
log T 



■ — rsi/Hi - +1 n 
- - [Si/H] = +2.0 
■ [Si/H] = +3.0 














••■ - — "* 




T = 9500 K 
eff 





-5 -4 -3 



T „= 15000 K 
eff 



-6 -5 -4 -3 -2 -1 
log T __ 



T „ = 1 1 000 K 

eff 




Fig. 1. Temperature differences between model atmospheres overabundant in Si by +1, +2 and +3 dex and the solar composition 
models for the set of effective temperatures. The differences are represented on the per cent scale AT/T with respect to the solar 
composition models. The optical depth scale logT ross is a logarithmic Rosseland scale. 



The combined influence of [Mn] = +3 and [Hg] = +5 on 
the model atmosphere structure does not differ from that of 
[Mn] = +3 alone. The maximum deviation is about 0.1% (or 
10K) for the highest T e s = 20 000 K, in the surface layers only. 
That is not surprising, since Hg belongs to the small changes 
group (Sect. 13X11 ). 

3.1.4. Special +2 dex case: Mg, Ca 

In this section we consider model atmospheres with chemical 
compositions which are not generally associated with CP stars 
(see Sect. 1 1.1 1 ). The reason for this is that, in fact, different chem- 
ical elements are of t en stratified in the atmo spheres of CP stars 
(Rvabchikova 2005; Rvabchi kova et alj|2005l) and the actual el- 
ement content in different atmospheric layers may be different, 
not necessarily the one deduced from the classic analysis (i.e. 
some layers may show non-typical enhanced or depleted values). 

At this point of our study, several chemical elements were 
assigned to the small changes group, according to the influence 
they have on the model atmosphere structure with the adopted 
abundance values. The fact is that some of these elements (Mg, 
Ca) are found to be stratified in the atmospheres of CP stars and 
can be strongly overabundant in some atmospheric layers. While 
it is clear that the elements not assigned to the small changes 
group (such as Si, Cr, Fe, etc.) are quite important in any varia- 
tions, the role of Mg and Ca if they are overabundant is unclear. 
So, we decided to consider these two elements overabundant by 
+2 dex (somewhat overestimated value, see notes about Ca and 
Mg in Sect. |3.2.Tll3.2.2| respectivelv) to complete the analysis of 
the small changes group elements. We consider these abundance 
values as a special case and treat them separately from typical 
underabundant values. They are marked with the asterisks (*) in 
Table |H 

Models with the excess of Ca can be assigned to the group 
of models with small changes in the temperature structure 
(Sect. 13.1.11 ) because only one model for T e ff = 11 000 K over- 
steps the 1% limit for the very surface layers (logr loss < -5.7) 
and reaches a 2% negative deviation. For high effective tem- 



peratures (T e g > 11 000 K), the models with [Ca] = +2 show a 
steady temperature increasing for the range of optical depths 
logT loss < -1 within the 0.5-1.0% bounds. 

The temperature changes for models overabundant in Mg 
fall into the group of the moderate effects (Sect. 13.1.31 ). For 
low effective temperatures (T e s < 11000K) the changes are 
still within the 1% threshold (except a narrow deviation by 
1.6% in between logr ross = — 1 and for T e s = 8000 K) while 
they exceed this limit for a wide range of optical depths with 
logr ross < -1 by about 2% (or 200 K) for the higher tempera- 
tures. 

3.1.5. Large (>3%) changes: Si, Cr, Fe and [M/H] 

In this section we consider the most interesting group of ele- 
ments (including scaled abundances) which produce the most 
noticeable temperature variations {AT/T > 3%) in the model 
temperature structure. 

Silicon. The comparison between temperature structures of the 
Si-peculiar and the solar composition model atmospheres are 
shown in Fig. Q] The main distinctive feature of the Si temper- 
ature profile is a broad plateau of increased temperature from 
about log T ross = upwards throughout the whole atmosphere 
for high values of the effective temperature (T e g > 11 000 K). 
For low effective temperatures (r er r < 11000K), there are no 
very remarkable changes in the model structure, with the excep- 
tion of the r er r = 8000 K, [Si] = +3 model, where the tempera- 
ture decreases considerably (even down to -4%) for all optical 
depths logT r0 ss < — 1. We discuss the reason for this behaviour 
in Sect. [3331 

The temperature deviations are about 0.5-1.5% (or 50- 
100K) for models overabundant by +1 dex, approximately 1— 
4% (or 100^00K) for +2 dex models, and about 1.5-6% (or 
100-800 K) for models with Si enhanced by +3 dex. As one can 
see in the figure, the "abundance - temperature deviation" rela- 
tion has a large spread in values due to its extreme dependance 
on the effective temperature. 
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Fig. 2. Temperature differences between model atmospheres overabundant in Cr by +1, +2 and +3dex and the solar composition 
models for the set of effective temperatures. The axes are the same as in Fig.Q] 



Chromium. Fig. |2] represents the temperature differences be- 
tween models with overabundant Cr and reference model at- 
mospheres. All subplots in this figure show similar temperature 
profiles with temperature increasing deep in the atmosphere and 
cooling in the upper layers. This is the usual temperature profile 
that model atmospheres with enhanced opacities are expected 
to demonstrate, due to the backwarming and surface cooling 
effects. All three abundance values produce progressively in- 
creasing changes in the atmosphere structure, cooling the optical 
depth regions with reduced temperature and heating the regions 
with increased temperature. 

The temperature deviation for the [Cr] = +1 model atmo- 
spheres is about 0.5% (or 50 K), for models with [Cr] = +2 its 
value is approximately 2% (or 100-200 K) and for model atmo- 
spheres peculiar in Cr by +3 dex the changes are about 4% (or 
200-500 K). The position (optical depth) of the inflection point 
separating the cooling and heating regions does not depend on 
the Cr abundance value and shifts evenly with the temperature 
increase from log r ross = -2 for T e g = 8000 K to log r ross = -3 
for r e ff = 20 000 K. Looking ahead, we can say that Cr seems to 
demonstrate the most organized temperature behaviour among 
all elements. 

Iron. The differences in the temperature structure between model 
atmospheres with overabundant Fe and the solar composition 
models are presented in Fig. [3] For comparison purposes, this 
figure also shows the temperature structure of the model atmo- 
spheres with scaled abundances. The influence of Fe is essen- 
tially the same as the influence of Cr: progressive heating and 
cooling of the deep and surface layers respectively as the Fe 
abundance increases. Generally, Fe demonstrates the most dra- 
matic changes among elements chosen to be considered in this 
section. We will discuss the Fe influence more in the following 
sections. 

The changes in the temperature structure for models over- 
abundant in Fe by +1 dex are about 1-3% (or 1 00-200 K) and 
for models with [Fe] = +2 those changes are about 2-6% (or 
200-800 K). 



Scaled abundances. The temperature structures of model atmo- 
spheres with scaled solar abundances are also shown in Fig. [3] 
As one can see these models produce much stronger temperature 
effects in comparison to the purely Fe-peculiar models. Only for 
the lowest effective temperature T e g - 8000 K are models with 
scaled abundances similar to those Fe-peculiar models. The fig- 
ure shows two evident heating and cooling regions in the model 
structure. The size of the heating region quickly grows, expand- 
ing towards the surface layers with increasing effective temper- 
ature. At the same time a weak cooling region shrinks more and 
more and even disappears for T e g = 20 000 K. 

Similar temperature patterns are sho wn by models with mag- 
netic line blanketing dKhan et alJ I2004L Fig. 1), and by mod- 
els with enhanced microturbulent velocity, most likely due to 
the fact that in all cases spectral lines of all chemical elements 
(rather than some of them) are strengthened. Another interest- 
ing detail is that [M/H] = +2.0 models provide almost exactly 
twice as much heating as [M/H] = +1.0 models do. Neither Fe 
nor Cr, Si or Mn show such a straightforward dependance, for 
any effective temperature. 

Model atmospheres with abundances scaled by +1 dex show 
changes in the temperature structure of about 2-5% (or 150- 
800 K), while for +2 dex models these values are around 3-10% 
(or 300-1600 K). 

3.1.6. Discussion on the temperature structure 

Here we want to summarize, explore and discuss in detail the 
main results of the temperature structure analysis and draw some 
conclusions for Sect. 13. 1 1 

After the analysis, within the range of abundances selected 
for consideration in this study (and which are typical for CP 
stars) we found that the most important elements to influence 
the temperature structure of model atmospheres of peculiar stars 
are Si, Cr and Fe. 

If we consider changes in the model atmosphere structure in 
the main line forming region (say from log r ross = -2 to 0) intro- 
duced by equally overabundant (+1 or +2 dex) elements of Si, 
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Cr and Fe taken separately, then, generally, the influence of Fe 
dominates over the influence of Si and Cr. The next most influ- 
ential element is Si, which in some cases (for high T e ff ) produces 
temperature changes which even slightly exceed those due to Fe 
(Teff = 9500 and 1 1 000 K, +1 dex) or are almost equal to them 
(e.g. r e ff = 13 000 K, +2 dex). Comparing Si-enhanced models 
to models overabundant in Cr by +3 dex shows that for low effec- 
tive temperatures, Cr produces a larger (T^tf - 8000 K) or simi- 
lar (T e ff = 9500 K) temperature deviation compared to Si, while 
for higher temperatures, Si overtakes Cr as was noted above. 
Consequently, putting iron, chromium and silicon in order of de- 
creasing temperature influence in the main line forming region 
we obtain the sequence of Fe-Si-Cr, where Si tends to reach or 
sometimes even exceed the level of the Fe influence for high ef- 
fective temperatures and sometimes yields to the influence of Cr 
for low effective temperatures. We find that this behaviour for Si 
is explained by its sharp temperature dependant influence pro- 
file noted above (Sect. 13. 1.5b . Having this element sequence in 
mind, later when we will consider their combinations in pairs 
we will test the influence of a "weaker" element on a "stronger" 
element (for example, how Si affects Fe-peculiar models or Cr 
affects Si-peculiar ones). 

An interesting result is that some models demonstrate an 
organized temperature behaviour with two distinct heating and 
cooling regions in comparison to the reference models. Elements 
that produce such temperature changes are Cr, Mn and Fe (which 
actually appear following one another in the periodic table of the 
elements). 

The heating region is located deep in the atmosphere and 
its magnitude clearly depends on the effective temperature 
and the element's overabundance value, steadily growing with 
the increasing of both of them for all mentioned elements. 
However, the temperature increasing pattern actually starts from 
Teff = 9500 K (not from the lowest temperature considered). 
This fact can be seen in either Fig. [2] for Cr or Fig. [3] for Fe, 
where the first subplots for = 8000 K show stronger heating 
than the second subplots for T e ff = 9500 K, while starting from 
the latter one, the temperature pattern seems to be held reliably. 
We did not provide a figure for Mn but it exhibits characteris- 



tics quite similar to those of Cr (including the dependence of the 
heating amplitude on the effective temperature). In addition to 
that, the scaled model atmospheres also demonstrate the same 
kind of behaviour, but even more pronounced (see Fig. [3]l in- 
cluding the exception for the T e $ = 8000 K value. 

The cooling region shows only the abundance dependant 
temperature behaviour. The higher the abundance is, the cooler 
the upper atmosphere. Besides that there are no other evident 
common traits for all elements to point out. 

Another feature for the group of Cr, Mn and Fe elements is 
that the inflection point which separates the cooling and heating 
regions progressively moves towards the surface with increas- 
ing effective temperature. For Mn and Cr, the position of this 
point does not depend on the abundance value but for Fe it does 
slightly and for high temperatures only {T e $ > 1 1 000 K). 

Model atmospheres with scaled solar abundances show tem- 
perature profiles that can not be simply replicated by overabun- 
dance of any single element separately (or otherwise). Even if 
some peculiar models are similar to those with the scaled abun- 
dances (e.g. models with [Fe/H] = +2.0 and [M/H] = +1.0 for 
T e ff = 9500 K in Fig. |3), there is no way to find the best model 
substitution in advance, without actually calculating one and 
looking for a fit, if any. We confirm that model atmospheres 
with scaled abundances can not be used to simply simulate the 
temperature effects of th e individual abundance patterns (see 
Pisku nov & Kupkal200U for a thorough discussion). 

We also tested the combined effect on the temperature struc- 
ture of the most influential elements. We calculated model atmo- 
spheres overabundant in two elements simultaneously (e.g. Fe 
and Cr, Fe and Si) considering all possible combinations among 
the values used in this study (see Table|2]i. 

Iron-chromium. We found that Cr enhanced up to +2 dex does 
not introduce much of a change in the temperature structure of 
the model atmospheres peculiar in Fe by +1 or +2 dex. Figure |4] 
illustrates this statement. The low sensitivity of the Fe-peculiar 
models to enhanced Cr content is explained by the fact that 
[Cr] < +2 dex itself does not change the atmosphere structure 
a lot - it produces less than a 2% effect (see Fig. |2j. What is 
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even more important is that model atmospheres peculiar in Cr 
have temperature profiles similar to the Fe-peculiar models as 
was discussed in Sect. 13.1.51 Altogether this results in a weak in- 
fluence of Cr enhanced by up to +2 dex on the Fe-peculiar mod- 
els, namely < 1 % (so the maximum influence is half of the Cr 
influence alone). 

As shown above, models with Cr overabundant by +3 dex 
show changes in the temperature structure which are quite large. 
When Cr overabundant that much is coupled with enhanced Fe 
there are obvious significant changes in the Fe-peculiar models 
due to Cr. However, it should not be forgotten that the [Fe] = +3 
value was not considered in this study since it is too high even for 
CP stars, thus comparison between models overabundant in Fe 
and in both Fe and Cr by +3 dex was not possible. We conclude 
that if Cr is less or equally abundant than Fe then its influence 
can be considered as insignificant (below 1%). 

Generally, model atmospheres with a larger Fe content are 
less sensitive to the excess of Cr. For instance, in the atmosphere 
with solar Fe abundance value the maximum temperature devi- 
ation resulting from the influence of Cr overabundant by +2 dex 
is 2% on average (see Sect. |3.1.51 >, while the atmosphere pecu- 
liar in Fe by +1 (+2) dex is changed by [Cr] = +2 by about 1.0 
(0.5)% respectively. In case a model is enhanced in Cr by +3 dex 
then the maximum deviation produced due to Cr on average is 
about 4% for the [Fe] = model, 2% for [Fe] = +1 and 1.5% for 
[Fe] = +2. 

Iron— silicon. Model atmospheres overabundant in both Fe and 
Si differ a lot from models enhanced only in Fe. The reason for 
this is the dramatically different effect of Si on the model atmo- 
sphere structure compared to that of Fe (compare Figs. Q] and 

Although the Fe-peculiar models additionally enhanced in Si 
by +1 dex are found to show less than 1% deviations with re- 
spect to those with no enhanced Si at all, the overall temperature 
changes which they show are very broad, covering almost the 
whole atmosphere. In case of [Si] = +2 or [Si] = +3 the con- 
trast between structures of the purely Fe-peculiar and FeSi mod- 
els is even more pronounced. Generally, an excess of Si seems 
to influence the temperature profile of the Fe-peculiar models 
in a unique way - uniformly increasing the temperature from 
log r ross = up to the surface, as it does for purely Si-peculiar 
models (see Sect. 13.1. BY 

The resulting effect of the exceptional influence of Si cou- 
pled with enhanced Fe is quite similar to what we see in the 
case of scaled solar abundances (compare the [Fe/H] and [M/H] 
profiles in Fig. [3). Indeed, our tests show that models peculiar 
in both Fe and Si by the same value (+1 or +2 dex) resemble 
the temperature profiles shape of the model atmospheres with 



scaled abundances ([M/H] = +1 or +2) fairly closely for low 
effective temperatures (T^ < 11 000 K). Of course, these FeSi 
models do not provide an exact fit but they are much closer 
than just Fe-peculiar models. If the Si abundance is increased 
more, by another +ldex, then the temperature profiles gener- 
ally provide better fit to those models with scaled abundances 
with high effective temperatures (T e ff > 1 1 000 K). For example, 
model atmospheres with [Fe] = +2 and [Si] = +3 exhibit a tem- 
perature structure that is very close (within < 1 % error bar) to 
the [M/H] = +2 model structure for r eff > 13 000 K. However, 
generally, despite this exceptional example, the discrepancy be- 
tween FeSi and scaled model atmospheres can be as much as 2— 
3% though it is still obvious that they are quite similar in shape. 

Finally, we consider how the changes introduced by Si de- 
pend on the Fe level, i.e. the sensitivity of the Fe-peculiar mod- 
els to additionally enhanced Si, in a similar way we did as for 
FeCr models in the previous subsection. Though, as it was men- 
tioned above the FeSi-peculiar models are substantially different 
from purely Fe-peculiar models, nevertheless there is a notice- 
able trend that the difference between them is reduced as the Fe 
content grows. Overall, this trend is 1.6-2.0 times less strong 
than that for Fe-FeCr models (i.e. an increasing amount of Fe 
does not decrease the Si influence as much as it does for Cr). 

Iron-silicon-chromium. The overabundance of Si, Fe and Cr si- 
multaneously by + 1 or +2 dex gives a better overall fit to models 
with scaled abundances than the FeSi-peculiar models do. But 
for high effective temperatures (T e ff > 1 1 000 K) additionally en- 
hanced Si is still needed, though maybe less in value. 

At this point we can conclude that the necessity of addition- 
ally overabundant Si may be attributed to the rest of the chemical 
elements whose influence grows with increasing effective tem- 
perature and works in a similar way as Si does. We find that 
among elements considered in this study, elements such as C, 
Mg and Ca would be good candidates to provide temperature 
increases in the upper atmosphere to fit temperature profiles of 
models with scaled abundances. 

It is interesting that for the lowest effective temperature 
(T e g = 8000 K), it is not found that Fe-peculiar model atmo- 
spheres additionally overabundant in Si or in both Si and Cr pro- 
vide a much better fit to models with scaled abundances than 
the purely Fe-peculiar models themselves. This happens pri- 
marily because the latter already closely resemble the temper- 
ature structure of the scaled abundance models (compare them 
in Fig. [3] for T e g = 8000 K and for the rest of T e ff values), and 
because enhanced Si and Cr result in small temperature varia- 
tions for this temperature. 
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Silicon-chromium. Model atmospheres peculiar in both Si and 
Cr probably have the most interesting element combination 
(among Si, Cr and Fe). The reason for this follows from the fact 
that Si influences the model structure in a different way than Cr 
does, and there is no substantially dominating element like Fe. 
Considering the structure of the SiCr models, it is clear that Si 
provides its distinct temperature profile which in turn is distorted 
by the influence of Cr. 

The main obvious effect produced by Cr is cooling of the 
upper atmosphere (see Fig. [2] for purely Cr-peculiar models, as 
an example). If Cr is overabundant by +2 dex then the tem- 
perature is decreased by about 1.5-2% in the surface layers 
(log r ross < -2) in comparison to Si-peculiar models, and the dif- 
ference reaches 2.5% in some cases. The Si-peculiar models en- 
hanced in Cr by +3 dex demonstrate about a 4% decreasing in 
the surface layers which reaches 5-6% sometimes. It is inter- 
esting that the cooling of the upper atmosphere due to Cr was 
not found to be dependent on the content of Si in any system- 
atic way. It stays almost the same in value, varying a bit in shape 
regardless of the Si value. 

The heating in the main line forming region provided by Cr 
to Si-peculiar model atmospheres is dependent on the effective 
temperature. For low effective temperatures (T e s < 1 1 000 K) its 
value does not depend on the Si content, most likely because 
Si itself does not change atmosphere structure notably for these 
temperatures (see Sect. 13.1.51 and Fig. [T). For higher tempera- 
tures there is a tendency that the Cr influence is being weakened 
with the increasing Si value. For example, the changes intro- 
duced by Cr overabundant by +3 dex in the model atmosphere 
with the solar chemical composition ([Si] = 0), on the average, 
are twice as large as those in the models with Si enhanced by 
+3 dex. 

3.2. Energy distribution 

In this section we study the energy distribution of the calcu- 
lated model atmospheres and continue the theoretical analysis 
of the depression at 520 A started in previous papers (see, e.g. 
iKhan & Shulvakl H)06a), which is frequently observed in the 
spectra of CP stars. 

Clearly, changes in the atmospheric temperature structure are 
directly related to the variations in the energy distribution. Now, 
after we have studied in detail the changes in the model atmo- 
sphere structure and found some regular behaviour in it, we do 
the same kind of analysis for the energy distribution. 

We are primarily interested in the analysis of those models 
whose temperature distributions show the largest deviation in 
comparison to reference model atmospheres to understand the 
interplay between them. However, we also analyse the energy 
distributions of models with weak deviations in the temperature 
structure. This analysis is performed in order to insure that small 
changes in the temperature structure are accompanied by small 
changes in the energy distribution and vice versa, i.e. to sup- 
port the initial element classification based on the temperature 
changes established in Sect. 13. II 

Generally, the changes in the energy distribution may be 
local (restricted to some wavelength range) or global (for the 
whole wavelength range, in practice the range adopted for model 
atmosphere calculation). The reasons for these changes in the 
case of the hydrogen-rich atmospheres (our case) are discussed 
next. 

Local changes appear as a result either of wavelength-limited 
spectral line absorption features (e.g. due to properties of the line 



density distribution of chemical elements which are important 
line opacity sources), or as a result of features in the contin- 
uum spectrum (discontinuities), which usually contribute essen- 
tially in a narrow UV region, and thus can be considered as local 
changes. These local changes are often responsible for global 
flux changes if located in the energetically effective regions (e.g. 
the far UV region for high T e ff) because of energy redistribu- 
tion: the excess or deficiency of energy in the local regions is 
balanced by a respective global flux deficiency or excess in the 
whole wavelength range. In other words, the global changes are 
the consequence of the local ones. 

Another possibility for global-like flux changes and tem- 
perature structure changes is variation of the mean molecular 
weight due to different chemical compositions. However in the 
hydrogen-rich atmospheres the actual influence of this variation 
is rather small and can be neglected. We will discuss this later 
with respect to the extreme cases of He-weak (Sect. 13.2.11 ) and 
Si-enhanced (Sect. [3~!2~3l model atmospheres. 

Before we consider flux changes in detail we note that, 
generally, global flux redistribution due to local features grows 
as the effective temperature grows, because the maximum of 
the energy distribution shifts towards ultraviolet, increasing the 
amount of energy being absorbed by the continuum features or 
by spectral lines in the far UV region. As a result there is usually 
a flux deficiency in the UV and a corresponding excess in the 
visual region. At the same time local flux features are tempera- 
ture dependant themselves, and the resulting energy redistribu- 
tion behaviour is a combination of both effects. 

We note also that hereafter if we describe any changes in 
the energy distribution without specifying explicitly its meaning 
we mean global changes (i.e. energy redistribution); the local 
changes we call flux features (e.g. the depression around 5200 A, 
a flux deficiency or discontinuity in the UV). 

The measurements of the global flux variations were per- 
formed primarily in the wavelengths interval 3000 through 
7000 A, i.e. in the region which is not largely congested with 
multiple absorption lines, unlike the far UV region, where the 
most local absorption occurs. 

The redistributed energy generally appears as a vertical shift 
with respect to the reference model atmosphere fluxes. However, 
this shift is not uniform over a wide wavelength interval, i.e. 
there is a continuum inclination with respect to the reference 
fluxes. This inclination is a key feature of the energy distribution 
of CP stars that makes it different (in the visual region it mimics 
fluxes produced by hotter and in the UV region by cooler normal 
stars). The question of the inclination of the continuum, which 
is closely connected to the calibration of photometric indices, is 
outside the scope of this work. It will be considered in detail as 
part of a forthcoming paper. 

For comparison purposes all fluxes were convolved with a 
gaussian with FWHM = 15 A, and comparisons were made on 
a logarithmic scale 2.5 • log F^ (i.e. magnitudes), where has 
units of erg s" 1 cm' 2 A -1 . 

Note that in the series of our papers about magnetic line blan- 
keting (see e.g. lKhan & S hulyak 2006a) we used units of log F^ 
to measure flux discrepancies. The label "magnitudes" thus has 
a different meaning in these earlier papers. 

Below, we use the same grouping of chemical elements that 
we used in Sect. 13.11 
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3.2.1. Small changes group 

The changes for model atmospheres peculiar in Sr, Eu, Hg, 
[Mg] = — 1, and [Ca] = -1 are less than 2-3mmag. For model 
atmospheres peculiar in He (for any considered values) the 
changes are about lOmmag. We found that the reason for these 
changes (both temperature and flux changes) for low effec- 
tive temperatures (r e ff < 1 1 000 K) is the variation (decrease by 
18.5%) of the mean molecular weight (recall that the tempera- 
ture changes due to depleted He are less than 0.3%). As we go to 
higher effective temperatures, both the mean molecular weight 
and an absorption edge of He i at 504 A are responsible for the 
changes in the temperature structure and flux with the growing 
influence of the latter one (recall again that He-depleted mod- 
els produce temperature changes less than 0.5% even for high 
temperatures, see Sect. 13.1 .TT ). All the models considered in this 
paragraph show very small changes in the temperature structure 
(which are well below 1 %, see Sect. 13.1.11 ), and so do the fluxes 
they produce. 

For model atmospheres with Ca overabundant by +2dex 
(special case model atmospheres which were assigned to the 
small changes group in Sect. 13-1 .4b , the flux changes are due to 
the flux discontinuity of Can at 1218 A and 1420 A. The max- 
imum local discrepancy between peculiar and reference model 
fluxes is about 4.0 mag deep, 100 A wide, located bluewards 
from the 1218 A feature, and appears for = 11 000 K (look 
at the temperature deviation from the reference model for this 
effective temperature, see Sect. 13.1.4b . The global changes for 
T e g = 11000K, however, are only about 15mmag in value. 
The maximum global flux variation of 25^-0 mmag appears 
for r eff = 15 000K, and for r eff = 20 000 K, it is about 10- 
25 mmag. 

Model atmospheres deficient in C and CNO demonstrate 
fluxes which differ from the reference model fluxes less than 
2-10 mmag for low effective temperatures (r e ff ^ 1 1 000 K) 
and reach values of 10-25 mmag for higher temperatures. The 
changes are due to flux discontinuities (absorption edges) in the 
UV region due to C i, at 1 101 A and 1240 A. 

3.2.2. Moderate changes group 

The changes in the flux distribution for model atmospheres pe- 
culiar in Ti are less than 3 mmag. We suppose that additional ab- 
sorption in spectral lines of Ti in the UV is responsible for this 
small flux redistribution. Although, in the temperature structure 
analysis (Sect. I3.1.3I I, Ti was assigned to the group of elements 
producing moderate changes, the weak variation in the energy 
distribution does not clearly support this assignment. The flux 
changes demonstrated by Ti -peculiar model atmospheres are ac- 
tually less than those produced by models peculiar in C or Ca 
(see the previous section). To explain this behaviour, we recall 
(see Sect. 13.1.31 1 that Ti-peculiar model atmospheres show the 
temperature deviation more than 1% only in the upper atmo- 
sphere (log t^s < _ 3) where there is little impact on the global 
flux distribution, since the main line forming region is well be- 
low this depth. Thus, Ti can now be assigned to the group of 
elements which provide small changes. 

Model atmospheres with enhanced Mn demonstrate energy 
distribution changes due to many absorption lines in the UV 
region. The redistribution effect grows with increasing effec- 
tive temperature. For models overabundant in Mn by +3dex, 
the changes in the energy are about 90- 120 mmag for high ef- 
fective temperatures (r e ff > 11 000 K) and are 50-70 mmag fol- 



lower temperatures. If we consider model atmospheres with 
[Mn] = +2 then these values are 40-50 mmag and 10-30 mmag 
for high and low temperatures respectively. 

Model atmospheres peculiar in Ni exhibit flux changes of 
about 25 mmag for low effective temperatures (r er j < 1 1 000 K) 
and about 30-60 mmag for higher temperatures. The changes are 
due to UV line absorption. 

Models overabundant in Mg by +2 dex (special case model 
atmospheres which were assigned to the moderate changes 
group in Sect. I3.1.4t show the strongest variation in the temper- 
ature structure among elements in the moderate changes group. 
This is a consequence of the influence of several continuum fea- 
tures of Mgi at 2515A and Mgn at 11 69 A and 1943 A. The 
global flux variation with respect to reference model fluxes is 
about 50-100 mmag on the average. The largest effect of about 
70-100 mmag occurs at T e g = 15 000 K. It is interesting that the 
model atmosphere with T s $ = 8000 K shows almost the same 
high values of the flux change because of a 1000 A wide (blue- 
wards from 25 15 A) and 2.0 mag deep continuum feature which, 
even for this low effective temperature, absorbs enough energy to 
produce the same effect as occurs in a model with almost twice 
the temperature. 

3.2.3. Large changes group 

Silicon. Model atmospheres with enhanced Si show substan- 
tial variation of the flux in their spectra. For [Si] = +1 these 
changes are about 10-40 mmag for low effective temperatures 
(T e g < 11000K) and about 25^40 mmag for higher tempera- 
ture. For [Si] = +2, these values are 60-170mmag and 80- 
170 mmag. Finally, for [Si] = +3 they are about 100-250 mmag 
and 200^-00 mmag, respectively. 

The reason for these changes is several large flux discontinu- 
ities due to Si i in the UV region (such as a feature at 1526 A due 
to ionization from the ground state, or a feature at 1682 A due 
to ionization from the first exited state) and strong absorption 
features in the UV region from ~ 1 100 to ~1600A due to Sin. 
We also found that there is a weak influence on the temperature 
structure and fluxes due to the variation of the mean molecular 
weight for Si-peculiar model atmospheres overabundant in Si by 
+3 dex. The content of Si (0.0295) is comparable with the con- 
tent of He (0.0483) for these model atmospheres and the mean 
molecular weight increases by +56% with respect to the solar 
composition model. However, we estimate that the resulting ef- 
fect on the temperature structure due to this variation is less than 
0.10-0.15% in the main line forming region and less than 0.6% 
for the upper atmosphere (log T loss < -4) in comparison to the 
reference models. In fluxes, molecular weight variations can be 
responsible for 15-30 mmag changes for the lowest and about 5— 
10 mmag changes for the highest effective temperature. Despite 
the fact that the increase of the mean molecular weight is rather 
large, the main effect on the model atmosphere structure appears 
as a variation of the atmospheric geometrical size - namely, the 
atmosphere contracts. However, as we mentioned above, there 
are only quite small changes in the optical depth scale. 

Here we recall the interesting feature found in the temper- 
ature distribution for the model atmosphere with r er r = 8000 K 
and peculiar in Si by +3 dex (see Fig. Q]). We find that the rea- 
son for this extraordinary temperature behaviour is that there are 
strong line absorptions and a number of distinctive flux discon- 
tinuities due to Si in the visual region (-4000-7000 A), which 
provide a unique broadband influence that is about 0.5 mag deep 
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with respect to the hydrogen continuum before the temperature 
correction procedure (in this case flux discontinuities give an 
example of non-local changes). These flux features are quickly 
weakened with increasing effective temperature and are almost 
absent for the next grid effective temperature value of 9500 K. 

There is also an important contribution due to Si spectral 
lines in the flux depression around 5200 A for low effective tem- 
peratures. For instance, for Si overabundant by +1 dex, there is 
a noticeable absorption contribution in the 5200-5500 A region 
for T e g = 8000 K. If Si enhanced by +2 dex then this contribu- 
tion is visible up to 1 1 000 K in the same wavelength range. For a 
value of +3 dex the line absorption in the 5200-5500 A region is 
clearly visible for effective temperatures up to T e $ = 13 000K. 
The reason for this feature is several strong Si i multiplets con- 
centrated near 5265, 5340, 5415 and 5495 A. 

Chromium. For Cr-peculiar model atmospheres, if Cr is over- 
abundant just by +1 dex, the flux changes are less than 10 mmag 
for low effective temperatures (T e g < 11 000 K), and less than 
25 mmag for higher effective temperatures. In case of [Cr] = +2 
the changes are 40-50 mmag and 40-70 mmag, respectively. 
Finally, if Cr is enhanced by +3 dex, then changes are 40- 
120 mmag for low effective temperatures and 100-150mmag for 
high ones. 

The energy redistribution happens due to many absorption 
lines of Cr in the ultraviolet region which effectively absorb en- 
ergy for the whole range of effective temperatures. 

The main contribution to the 5200 A depression is because 
of the absorption in the region of 5 150-5450 A due to numer- 
ous spectral lines of Cr. The effect grows with the growing Cr 
abundance and decreases with increasing effective temperature. 
The contribution can be seen clearly up to T e g = 11 000 K if the 
highest overabundance value of [Cr] = +3 is adopted. 

Iron. Model atmospheres peculiar in Fe by +ldex show flux 
changes of about 50 mmag for low effective temperatures 
(Teff < 1 1 000 K) and 70-120 mmag for higher effective temper- 
atures. For models with Fe overabundant by +2 dex, these values 
are 120-170mmag and 200-370mmag respectively. 

There are several flux discontinuities in the UV region due 
to Fe, but they only appear for the lowest effective temperature, 
and thus do not contribute substantially to the energy redistri- 
bution. The global flux changes are primarily because of energy 
absorption in the UV region by a great number of Fe lines for 
all effective temperatures. We found that the influence of the al- 
tered value of the mean molecular weight (+11% for [Fe] = +2) 
is absolutely negligible. 

Also, Fe contributes significantly to the 5200 A depression 
due to a great number of Fe spectral lines absorbing energy in a 
wide wavelength range from 5000 to 5500 A. This contribution 
has a weak tendency to be stronger a bit bluewards of 5200 A. 
If Fe is enhanced by +1 dex then the depression is seen for ef- 
fective temperatures up to 15 000 K. If Fe is overabundant by 
+2 dex then the depression is clearly recognizable up to the high- 
est effective temperature of 20 000 K, gradually decreasing with 
increasing effective temperature. 

Scaled abundances. Model atmospheres with a chemical com- 
position scaled by + 1 dex produce fluxes which differ from 
the reference model atmosphere fluxes by 70-200 mmag for 
low effective temperatures (T e ff < 11 000 K) and by about 120- 
250 mmag for higher temperatures. For model atmospheres with 
[M/H] = +2, flux changes are about 120-400 mmag for low ef- 
fective temperatures (T e g < 1 1 000 K) and about 350-500 mmag 
for higher effective temperatures. 



Concerning the flux depression around 5200 A, there is 
clearly a combined contribution due to Si, Cr and Fe to this de- 
pression. Putting together their influences as outlined above, we 
conclude that Fe is the principal contributor for the whole range 
of effective temperatures, while Cr and Si are important primar- 
ily for low effective temperatures. 

3.3. Photometry (the uvby/3 and Aa systems) 

Photometry provides us with a powerful and well-developed tool 
for a preliminary analysis of observed stellar fluxes, capable of 
rapid measurement of a large number of stars. Quite often the 
photometric indices are the only available data about the stel- 
lar energy distribution, which makes the theoretical photometry 
analysis especially valuable. 

Various photometric systems have been designed to work in 
different spectral intervals to probe specific broad- or narrow- 
band spectral features. In this section we consider two photo- 
metric systems. The first one is the Stromgren-Crawford uvby/3 
photometric system which is a general one. The second system 
is the peculiarity index a designed to measure the strength of the 
5200 A depression with respect to normal stars (the Aa photo- 
metric system). 

We computed synthetic photometric colors in these two sys- 
tems following the proc edure used in our previous studies (e.g. 
Khan & Shulvak 2006a b). There is also a detailed description 
of how to calculate the synthetic a index by Kup ka et al. I J20031 
Sect. 2.1). 

An in-depth analysis of the photometric indices is beyond the 
scope of this paper, thus here we restrict ourselves to a summary 
of the photometric changes found. First we present a general dis- 
cussion of the photometric index values for various models, than 
we describe the a vs. b - y diagram specifically for model atmo- 
spheres of the large changes group. All values were rounded to 
1 mmag accuracy. 

A thorough investigation of the spectral energy distribution, 
photometric colors and calibration of photometric systems will 
be presented in a forthcoming paper. 

3.3.1. Discussion on the photometric indices 

The results in this section are presented with respect to the ref- 
erence model atmospheres. 

Small changes group. On the whole, the variations of all pho- 
tometric indices of the uvby/3 system in this group are approxi- 
mately within a ±3 mmag limit. However there are quite a num- 
ber of exceptions to this rule which can be rather large, primarily 
for ;«i for the lowest effective temperature (8000 K) and for c\ 
for all temperatures. The reasons for this are various line ab- 
sorption features (such as H and K Ca 11 lines or Sr 11 and Eu 11 
lines) and the global energy redistribution. An important thing 
to notice is that none of the models produce changes in b - y of 
more than +2 mmag, except those with [Ti] = +2 for low effec- 
tive temperatures (e.g. A(b - y) — +13 mmag for T e ff = 8000 K). 
Finally, no models show variations in the peculiarity index a of 
more than +2 mmag. 

Moderate changes group. It is hard to point to any particular 
local features contributing to a certain photometric filter which 
cause changes to the derived photometric indices for models in 
this group. We suppose that the main cause of the variations is 
the global energy redistribution which affects values obtained in 
all photometric filters. The indices of the uvby/3 system show 
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positive as well as negative variations, while the most extreme 
values are for low effective temperatures. The temperature in- 
dicator b — y exhibits deviations of about +4 mmag with some 
exceptions (e.g. mean deviation of about -15 mmag for models 
with [Mn] = +3). The variations for a are about ±5 mmag show- 
ing extreme values of + 15 mmag for the lowest r e ff. 

Large changes group (Si, Fe, Cr and [M/H]). The global energy 
redistribution and local flux depressions are especially strong for 
this group of elements and are responsible for changes in all in- 
dices. 

All the models show positive changes in the m\ index 
with respect to the reference model atmospheres, with some 
exceptions (up to -18 mmag) for purely Cr-peculiar models 
([Cr] < +3). The model atmospheres which are peculiar in one 
or several elements by no more than + 1 dex show typically about 
5-20 mmag variations which can rise to 20-50 mmag (primarily 
for r e ff = 8000 K). If one or several elements are overabundant 
by +2 dex then these changes are usually about 10-100mmag 
(with largest values of 30-170mmag). Finally, if one or several 
elements are enhanced by +3 dex then changes are about 20- 
120mmag, with maximum values of 40-180mmag. Note that 
the smallest Ami values are for the purely Cr-peculiar models. 

Generally, all the models demonstrate negative variations 
in the ci index (i.e. the Balmer jump becomes smaller), 
but there are some exceptions for low effective temperatures 
(T e ff < 9500 K) showing positive values for Si- and Cr-peculiar 
models (up to +32 mmag). The model atmospheres which are 
peculiar in one or several elements by no more than + 1 dex 
show about -5 to -50 mmag variations, which peak to -10 
to -70 mmag. If one or several elements are overabundant by 
+2 dex then these changes are about -15 to -200 mmag (with 
peaks at -20 to -300 mmag). Finally, if one or several ele- 
ments are enhanced by +3 dex then changes are about -30 to 
-200 mmag with peaks at -40 to -450 mmag. Note that the 
highest Aci values are for the purely Cr-peculiar models. 

For the (3 index, positive A/3 values tend to occur at the low 
end of the effective temperature scale while negative values oc- 
cur at its high end. The majority of the values are concentrated 
within a +10 mmag range and peak up to about ±25 mmag for 
models strongly peculiar in Si for low effective temperatures. 

The b — y and a indices are considered in the next section. 

3.3.2. The diagram of the peculiarity index a vs. b - y 

For the second part of the photometry analysis we consider the 
diagram of a vs. b - y. Thus we can explore the Aa photometric 
system the way it was designed, i.e. with respect to the normal 
stars of the same temperature indicator color (the b—y index in 
our case). 

Figure |5]represents model atmospheres of the large changes 
group computed in this work (except those with strong overabun- 
dance values of [Si] = +3 and [Fe] = +2, to decrease number 
of points in the plot and provide better picture resolution and 
thus visualization). Note that the conclusions drawn below (il- 
lustrated on the basis of the limited sample in Fig.|5]l are for the 
whole parameter space. We advise the reader to read thoroughly 
the caption of the figure for details. The analysis revealed several 
interesting features: 

1 . An important feature to notice is that none of the models pe- 
culiar in Fe (the plus (+) symbols) show negative Aa values, 
they all appear above the normality line. 

2. Model atmospheres which show negative Aa values have 
T e g > 1 1 000 K and are primarily concentrated leftwards 
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Fig. 5. Diagram of the peculiarity index a vs. the effective 
temperature indicator b—y. The figure represents model atmo- 
spheres of the large changes group (except models with strong 
abundances of [Si] = +3 and [Fe] = +2). The non-peculiar mod- 
els are model atmospheres assigned to the small changes group 
(Sect. [XTTTT ). and indicated on the plot with dots (•). All models 
which are peculiar in Fe by +1 dex are shown with a plus sign 
(+), i.e. this includes purely Fe-peculiar, scaled, FeCr-, FeSi- 
and FeSiCr-peculiar model atmospheres. Models peculiar in Cr 
and Si only are shown with the square (□) and the asterisk (*) 
symbols, respectively. The upward pointing triangles (A) repre- 
sent SiCr-peculiar model atmospheres. The black solid line is 
the normality line arj of the peculiarity index Aa - a - ao, and 
is a parabolic least squares fit through points which correspond 
to the reference model atmospheres. The inclined dotted lines 
represent isolines of the Aa index spaced 5 mmag. 



from b — y = -0.045 (the value shown by the reference 
model atmosphere of T e ff = 13 000 K). All of these models 
have the solar content of Fe. 
3. The locus of points in the right upper section shows the regu- 
lar highly organized pattern which we investigate below. The 
same feature is also clearly recognizable in the middle of the 
figure leftwards from the vertical zero line. 

To understand the regular pattern found in Fig. [5] we plotted 
three separate b - y vs. a diagrams for models of different effec- 
tive temperatures (Figs. [6HH}- We have chosen values of 8000 K, 
9500 K and 15 000 K because the reference model atmospheres 
of these effective temperatures on the b—y scale represent the 
left, middle and the right parts of the diagram shown in Fig. [5] 
Thus, this way we investigate two patterns found above (point 
number 3) and investigate whether there is the same regular be- 
haviour for the high temperature models which is harder to iden- 
tify within the diagram in Fig. [5] 

Figure [6] shows the b — y vs. a diagram for T e g = 8000 K. 
The caption for the figure provides the explicit description of 
symbols and labels. The figure demonstrates that there are three 
main directions (or axes) on the diagram associated with grow- 
ing content of each of the elements Si, Cr and Fe. Moreover, 
there is an obvious hint that all the models tend to be located in 
the intersections of the isolines. For example, if we sequentially 
follow along the Si, Cr and Fe axes by the amount correspond- 
ing to a + 1 dex change in each move then we get to the point 
where the model atmosphere overabundant in Si, Cr and Fe by 
+ 1 dex is located (the one marked with the symbol ffi in the right 
bottom part of the diagram). Consider another example, of the 
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Fig. 6. Diagram of the peculiarity index a vs. b-y for all 
model atmospheres of the large and small changes groups with 
Teff = 8000 K. Models are represented with the same symbols as 
in Fig. [5] The plus sign (+) retained for all Fe-peculiar models, 
while the crossed circle (ffi) is for SiCrFe models and the star 
symbol (★) is for models with scaled abundances. Thick lines 
with numerical labels in dex units represent peculiarity axes of 
Si, Cr and Fe, e.g. axes formed by model atmospheres peculiar 
only in one of these elements. For visualization purposes there 
are isolines on the plot: dashed lines are the isolines of the Fe- 
axis, the rest of isolines are the solid lines. All curves are least 
squares parabolic or linear fits through the corresponding points. 



model atmosphere overabundant in Si by +2 dex and in Cr by 
+3 dex. Again, to get to the point on the diagram which repre- 
sents this model we should follow the Si direction to the point 
labelled +2 dex and then along the Cr direction to the point with 
the +3 dex mark. 

Of course, there are not perfectly concordantly aligned iso- 
lines everywhere on the diagram. For instance, in the left part of 
the diagram, there is a pattern break for model atmospheres pe- 
culiar in Si by +3 dex. The distance between the Si-axis marks 
stretches slower than the corresponding one between nodes on 
the isolines of the SiCr-peculiar models (see connected triangles 
(A) symbols along the Si direction on the diagram). 

Similar diagrams for T e $ = 9500 and 15 000 K are shown in 
Figs. [7J and [8] respectively. The analysis of all these figures and 
those for other effective temperatures reveals some interesting 
properties: 

1 . As the effective temperature grows the low curvature Si-axis 
gradually rotates counter-clockwise sinking below the nor- 
mality line at around 10000K. 

2. The Cr-axis shows the same behavior, going below the nor- 
mality line at around 13 000 K. At the same time it changes 
its shape from the initially straight line to the arched one. 

3. Finally, in contrast to items above, the Fe-axis is very stable 
in its shape and inclination. In fact, it is almost straight and 
virtually perpendicular (on the scale of Fig. |6]l to the nor- 
mality line for any effective temperature. If all the Fe axes 
for different effective temperatures are plotted on the same 
diagram then they appear to be almost perfectly parallel. 

This analysis supports our earlier conclusion that Fe is the 
main contributor into the 5200 A depression for the whole range 
of effective temperatures (Sect. [3~!23T >. 
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Fig. 7. The same diagram as in Fig. |6]but for = 9500 K. Dash 
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triangles A) may not be a fully adequate approximation in this 
particular case. 



0.645 








0.635 


\* Xv 


: T eff = 


15000 K 


0.625 








djO.615 

CO 

E 














to 0.605 








0.595 








0.585 








0.575 






/ +3 
/Si 



-0.16 



-0.14 



-0.12 -0.1 
b-y, mag 



-0.08 



-0.06 



Fig. 8. The same diagram as in Fig. [6] but for T e ff = 15 000 K. 
Dash lines in this plot are the isolines of the Si-Fe grid. 



It also demonstrates that Fe is indeed the key element for the 
fundamental property of the Aa system to show positive values 
for CP stars with overabundant Fe-peak elements in their atmo- 
spheres, and negative values for CP stars with underabundances 
of the Fe-peak ele ments (A Bootis star s) (see the summary on Aa 
measurements bv lPaunzen et al.ll2005l Table 3). 

The diagram analysis confirmed the role of Cr and Si as im- 
portant contributors into the magnitude of the 5200 A depres- 
sion (measured through the Aa system) for low effective temper- 
atures. 

To extend our study and ensure the reliability of the diagram 
properties outlined above we computed several additional model 
atmospheres for the effective temperature of 7000 K and found 
indeed for them the same features. 

Also we found that model atmospheres scaled by + 1 dex 
show b-y values which differ from those of the reference model 
atmospheres by a small to modest amounts: -1, -19, -24, -25, 
—21, -lOmmag (in order of increasing effective temperature). 
These val ues are actually in a perfect agreement with numerical 
results bv lKuruc3 (jT993a), and their modest values partly explain 
why we did not find any systematic errors in the determination 
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of the fundamental parameters of magnetic CP stars with scaled 
abundances in their atmospheres (Khan & Shulyak 2006a). 

At the same time Figs. |6]-[8] show that if the model atmo- 
sphere is overabundant in Si, Cr or Fe (separately or together) by 
more than + 1 dex then this may result in a quite large deviation 
of the temperature indicator b - y, that is explainable by the flux 
redistribution effects (see Sect. 13.2.3b . However, in some cases 
simultaneous overabundance of two or even three of these ele- 
ments produces less variation of the b — y value than they would 
separately. 

3.4. Hydrogen line profiles 

To study the effects due to abundances other than solar on the 
hydrogen line profiles we calculated profiles of the Ha, H/3 and 
Hy lines and compared them with those for reference model at- 
mospheres. 

In accordance with the results obtained in Sect. 13.11 (model 
atmosphere temperature structure) we decided that model at- 
mospheres with peculiar abundances such as [M/H] = +1, 
[Fe] = +1 and [Si] = +2 would be sufficiently representative for 
the analysis of the hydrogen line profiles. The comparison be- 
tween models with an effective temperature of 1 3 000 K, which 
we find to be the best case to illustrate the following discussion, 
is represented in Fig. [9] 

We found that hydrogen line profiles for models with 
[M/H] = +1 and [Fe] = +1 are essentially the same for high 
effective temperatures (T^ > 11 000 K) and show a noticeable 
negative deviation from the reference line profiles in the re- 
gion about +5 A around the line center. Moreover, the differ- 
ence between peculiar and reference line profiles almost does 
not change as the effective temperature changes. For low effec- 
tive temperatures (T^ < 9500 K), the line profiles for the pe- 
culiar models do show broad negative variations (especially for 
8000 K) from the reference line profiles and differ from each 
other by about 1-2% (models with scaled abundances produce 
a larger effect). Overall we find that the maximum variation in 
profiles of Ha, H/3 or Hy is about -4% (generally, the strongest 
variation is for Ha, the weakest for Hy). 

Model atmospheres peculiar in Si show primarily nega- 
tive and broad line profile variations for low effective temper- 
atures. These changes are about -1% to -2% for models with 
r eff = 8000 K and about -3% for T eff = 9500 K. For higher ef- 
fective temperatures the differences between line profiles of pe- 
culiar and reference model atmospheres are still broad but they 
become positive showing variations of up to — 1-2%. 

3.5. Abundance analysis 

In this section we present results of a brief investigation on abun- 
dance analysis in the atmospheres of CP stars. An introduction 
to the problem is given in Sect. 11.21 We study the effects of the 
individual abundance patterns on the abundance determination 
procedure through the spectrum synthesis method. The idea of 
the method is to fit theoretical calculations to the observations. 
Usually, this involves the initial choice of the model atmosphere 
(generally calculated for a scaled solar chemical composition) 
which is believed to be the best approximation, and the follow- 
ing use of a spectrum synthesis code to compute the high res- 
olution spectrum. The atmospheric abundances are to be varied 
during the synthesis procedure to match the synthetic spectrum 
to the observed one. 



The questionable point of this method is that the model at- 
mosphere which is used for the spectrum synthesis is considered 
as a parameter, not as a function. In other words, the model is 
not recalculated in accordance with the results of the abundance 
analysis. 

We study the effects of this approach as applied to the abun- 
dance analysis of Si, Cr and Fe (elements of the large changes 
group, Sect. 13.1. BY The study covers abundance values of Si and 
Fe up to + 1 dex and Cr up to +2 dex. 

First, we calculate a reference spectrum for a particular 
model atmosphere from our grid (Si-, Cr- and/or Fe-peculiar 
model). Second, we use model atmospheres for scaled chem- 
ical compositions ([M/H] — and [M/H] = +1) to synthesize 
two fitting spectra, assuming the same abundances of Si, Cr or 
Fe that were used to compute the reference spectrum. Then, we 
compare each of the fitting spectra to the reference one and, if 
necessary, vary the abundance values of Si, Cr or Fe in order to 
fit the reference spectrum. The changes in abundances required 
to fit the reference spectrum reflect errors or other unwanted, 
possibly systematic, effects following from the use of the incor- 
rect model atmosphere. All deviations are measured with respect 
to the reference spectrum. 

The spectral line data for the spectrum synthesis procedure 
was extracted from VALD. We used the same lines of Si, Cr 
and Fe which were adopted for t he s tratification analysis of 
HD 204411 by iRvabchikova et al.1 (l2005l Table 2). The advan- 
tage of this compilation is that all the lines are carefully se- 
lected, representing samples of the low- and high-excitation lines 
of neutral and once ionized atoms. The data cover the spectral 
range from -5000 A to ~6400 A and consists of 23 lines of Cr, 
23 lines of Fe and 9 lines of Si. All spectra were computed with 
zero microturbulence and no rotation. 

Here we note that in our analysis we suppose that the hy- 
pothetical analyst knows precisely the fundamental parameters 
of the stars being analyzed, although this is not correct for real 
stars, and especially not for CP stars. There are usually also other 
uncertainties and difficulties in working with real spectra which 
influence the actual abundance analysis procedure. Our theoret- 
ical study is not intended to simulate accurately the problem the 
observer faces, but to answer the question of what are the am- 
biguities of the abundance analysis (including the analysis of 
vertical stratification) using classic model atmospheres, if this 
analysis is performed in the ideal conditions of other variables 
being fixed and known. For the same reason of providing ho- 
mogeneity to the study, we use the same spectral lines for all 
effective temperatures. 

3.5.1. Fe abundance 

We use the Fe-peculiar model atmosphere with [Fe] = +1 to 
calculate the reference spectrum. The comparison of the fit- 
ting spectrum ([M/H] = 0) with the reference one showed 
monochromatic flux discrepancies of about +2% on average, 
with positive (+4%) and negative (-3%) peaks. Despite these 
small values we were able to find some clear systematic effects, 
which actually appear to be common properties for other spectra 
and other elements tested within this abundance analysis. Thus 
we describe these effects and give some examples as they apply 
to this particular fit: 

1. The spectral lines of neutral and once ionized atoms show 
different deviations with respect to the reference spectrum, 
requiring different abundance corrections to fit them. 
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Fig. 9. Comparison between Ha, H/? and Hy line profiles computed for model atmospheres with T e ff = 13 000K and different 
chemical composition (peculiar abundances [M/H] = +1, [Fe] = +1, [Si] = +2 and the reference solar abundances [M/H] = 0). 
The upper subplots show the difference (A/ = / - 7[m/h]=o) between line profiles calculated for models with peculiar abundances 
and those for the reference model atmospheres. The normalized line profiles {I/I c ) are represented in the bottom subplots. 



For example, for T e g = 11 000 K, the Fe i lines show negative 
while the Fe n lines positive variations that require negative 
and positive corrections to the Fe abundance, respectively. 

2. It is impossible to fit the core and the wings of strong lines 
simultaneously with the same abundance value. This effect 
decreases as the effective temperature grows (and the number 
of strong lines goes down). 

For example, for T e ff = 8000 K, all the Fei lines show this 
effect, demonstrating more than the +3% deviations in the 
core and wings. 

3. The low- and high-excitation lines and the lines of different 
intensities require different abundance corrections. 

For example, for T e g = 15 000K, the high-excitation Fen 
lines require an abundance value which is ~0. 1 dex higher 
than that for low-excitation lines. 

The comparison of another fitting spectrum ([M/H] = +1) 
to the reference one demonstrates quite similar behaviour except 
that all changes are of the opposite sign. 

Since we do not expect that only Fe would be overabun- 
dant in the real atmosphere of a CP star, we calculated a sec- 
ond reference spectrum using a model atmosphere peculiar in 
Si, Cr and Fe by +1 dex. The comparison of the fitting spectrum 
([M/H] = 0) to the new reference one shows +3.5% monochro- 
matic flux variations on average with positive (+5%) and nega- 
tive (-6%) peaks along with the same properties found above. 

Another fitting spectrum ([M/H] = + 1) shows an almost per- 
fect agreement (about 0.5%) with the new reference spectrum in 
the range of effective temperatures from 9500 K to 13 000 K due 
to a good agreement between the temperature structures of these 
two models (see Sect. 13.1. 6l . For other T e g values the deviation 
is about 1-2%. 

Finally, we find that the ambiguity in the determination of the 
Fe abundance due to use of various lines is within a ±0.25 dex 
error bar, and apparently depends on the differences in the at- 
mosphere structure (i.e. density, pressure, temperature and other 
quantities influencing line profiles) of the model being tested. 



3.5.2. Cr abundance 

As in the previous subsection we adopted two models peculiar 
in Cr to calculate the reference spectrum. One model is peculiar 
only in Cr by +2 dex, another one is more realistic and addi- 
tionally peculiar in Fe by + 1 dex. The latter appears to produce 
the best agreement in the sample (about 1-2% discrepancies) if 
compared to the [M/H] = +1 fitting spectrum. 

During the analysis we revealed all the same systematic ef- 
fects outlined for Fe. The uncertainties found in the analysis are 
also about the same values as for Fe, i.e. inside a +0.25 dex error 
bar. 



3.5.3. Si abundance 

Again, we used two model atmospheres to compute the refer- 
ence spectrum. The first one is the Si-peculiar model ([Si] = + 1), 
the second one is the SiCrFe-peculiar model with all these el- 
ements overabundant by + 1 dex, i.e. the same model we used 
for the abundance analysis of Fe. The best agreement (about 
1% deviation) is found for the [M/H] = fitting spectrum and 
the reference spectrum of the Si-peculiar model, and for the 
[M/H] = +1 fitting spectrum and the reference spectrum of the 
SiCrFe-peculiar model. Overall, the same systematic effects and 
approximately the same error bar were found for Si as for Fe and 
Cr. 



3.5.4. Summary on the abundance analysis 

In spite of the fact that the typical error bar found for abundance 
analysis in the visual region is about +0.25 dex, the particular 
error value depends on the model atmosphere being used for the 
analysis, and varies for different effective temperatures. In fact, 
some models provide a very good fit to the reference spectrum 
with almost no error, while others result in differences beyond 
this error bar. Apparently, these discrepancies follow from the 
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differences in the temperature structure of the model used for 
analysis and the actual atmosphere (see Sect. 13.1.61 . 

Also, considering the results of our numerical experiments 
as applied to stratification an alysis using homogeneo us model 
atmospheres (see for instance Rvabchi kova et al.ll2005l) . we find 
that the uncertainty of the value of the vertical abundance gradi- 
ent is within an 0.4 dex error bar. 

4. Conclusions 

We have computed and analyzed a grid of model atmospheres 
of chemically peculiar (CP) stars. The main purpose was to per- 
form a systematic homogeneous analysis of the effects of the 
individual abundance patterns on the model atmosphere struc- 
ture, energy distribution, photometric indices (in the uvby/3 and 
Aa systems), hydrogen line profiles, and on an abundance deter- 
mination procedure. 

The grid of model atmospheres consists of more than 300 
models and represents the following types of stars: CP1 (Am 
stars), CP2 (Si, Cr-Sr-Eu A and B stars), CP3 (Hg-Mn, B-type 
stars) and CP4 (He-weak, B-type stars). All models were cal- 
culated with the LLmodels code (version 8.4) assuming no mag- 
netic field and no convection to highlight abundance effects only. 
Consequently, we considered model atmospheres with different 
chemical composition and compare them to the reference mod- 
els with solar abundances. Detailed data about the chemical el- 
ements examined and the adopted abundance values are repre- 
sented in Table [2] 

The conclusions and main results of the study are summa- 
rized below: 

1 . The majority of the tested chemical elements (within the lim- 
its of abundance values considered) produce less than 1% 
variations in the model atmosphere temperature profile and 
fall into the small changes group. These elements are He, 
C, N, O, Mg (deficient), Ca, Sr, Eu, Hg. According to the 
results, a 1 % limit was adopted as the error bar threshold. 

2. Several elements were assigned to the moderate changes 
group in accordance with the effects they produce on the 
model atmosphere structure (1-3%) and the energy distri- 
bution. These elements are Mn, Ni, Mg (enhanced). 

3. The group of elements which produce large changes in the 
model atmosphere structure (more than 3%) and energy dis- 
tribution consists of Si, Cr, Fe and scaled abundance patterns. 

4. If we consider changes in the temperature structure produced 
by the elements Si, Cr and Fe in the main line forming region 
then we find that the sequence of Fe-Si-Cr represents these 
elements from the most to the least influential. 

5. Model atmospheres peculiar in Cr, Mn and Fe demonstrate a 
very similar, highly organized, temperature behaviour. There 
are two distinct cooling and heating regions in the upper and 
the lower (i.e. the main line forming region) atmosphere, re- 
spectively. As the abundance value grows, the temperature 
drops in the upper atmosphere and increases in the lower 
atmosphere. The inflection point which separates these two 
regions moves outwards with growing T e ff. The magnitude 
of the heating region steadily grows with increasing T e s ex- 
cluding the lowest value of 8000 K. The same effects, except 
for the cooling region feature, apply to models with scaled 
abundances as well. 

The reason why Cr, Mn and Fe demonstrate such behaviour 
is likely the similarity of their energy level configurations 
(and the corresponding spectral line distribution patterns), 
and their relatively high content in the atmosphere. 



6. We concluded that the elements Si and Fe are the main 
providers of two different types of the temperature changes 
to produce the same distinctive temperature structure that 
model atmospheres with scaled abundances do. In other 
words, the temperature structure of a model atmospheres pe- 
culiar only in Si and Fe may be close to the temperature pro- 
file of a models with scaled abundances. The agreement is 
quite good for low effective temperatures, while for high ef- 
fective temperatures (T e s > 11 000 K) the cumulative effect 
of all other chemical elements is required. We suspect that 
overabundant C, Mg and Ca are of the most significance. 

7. We confirm that model atmospheres with scaled abundance 
patterns can not be used to simulate accurately effects of the 
individual abundance patterns. 

8. We found that Fe is the principal contributor into the 5200 A 
depression for the whole range of effective temperatures, 
while Cr and Si are important primarily for low effective 
temperatures. 

9. The analysis of the diagram of the peculiarity index a vs. 
b — y for the model atmospheres peculiar in Si, Cr and Fe 
revealed regular patterns in the locus of points represent- 
ing those models. In fact, there are three separate directions 
(axes) associated with the growing content of each element 
on the diagram. The inclination of the Fe-axis to the nor- 
mality line flo and the fact that it does not depend on the ef- 
fective temperature clearly demonstrate that Fe is indeed the 
key element to the fundamental property of the Aa system to 
recognize CP stars with overabundant (Aa > 0) and under- 
abundant (Aa < 0, i.e. for A Bootis stars) Fe-peak elements 
in their atmospheres. 

10. We investigated an abundance analysis procedure based on 
theoretical atmospheres with individual abundance pattern 
using models with scaled solar composition and spectral 
lines in the visual region. We find that the error bar of the 
analysis, which occurs in the context of the ideal conditions 
of a theoretical study, is of order 0.25 dex, and that the par- 
ticular error value depends on the model atmosphere being 
used for such analysis, and varies for different effective tem- 
peratures. 

1 1 . Considering the results of our numerical experiments as they 
apply to the stratification analysis using homogeneous model 
atmospheres, we conclude that uncertainty of the value of the 
vertical abundance gradient is within an 0.4 dex error bar. 
That, of course, does not eliminate the phenomena of the 
abundance stratification in the atmospheres of CP stars, how- 
ever it can effectively increase or decrease the gradient value 
deduced from the observations. That, in turn, may be essen- 
tial in comparison with the results of self-consistent diffusion 
calculations. 

In a forthcoming paper, we are planning to focus on the study 
of observed energy distributions of CP stars, their photometric 
indices, and the calibration procedure for the determination of 
fundamental parameters. 
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